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Abstract. 

Inspired by a previous work by McClure and Dyer (Class. Quantum Grav. 23, 1971 
(2006)), we analyze some solutions of the Einstein-Maxwell equations which were originally 
written to describe charged black holes in cosmological backgrounds. A detailed analysis 
of the electromagnetic sources for a sufficiently general metric is performed, and then 
we focus on deriving the electromagnetic four-current as well as the conserved electric 
charge of each metric. The charged McVittie solution is revisited and a brief study of 
its causal structure is performed, showing that it may represent a charged black hole in 
an expanding universe, with the black hole horizon being formed at infinite late times. 
Charged versions of solutions originally put forward by Vaidya (Vd) and Sultana and 
Dyer (SD) are also analyzed. It is shown that the charged Sultana-Dyer metric requires a 
global electric current, besides a central (spherically symmetric) electric charge. With the 
aim of comparing to the charged McVittie metric, new charged solutions of Vd and SD 
types are considered. In these cases, the original mass and charge parameters are replaced 
by particular functions of the cosmological time. In the new generalized charged Vaidya 
metric the black hole horizon never forms, whereas in the new generalized Sultana-Dyer 
case both the Cauchy and the black hole horizons develop at infinite late times. A charged 
version of the Thakurta metric is also studied here. It is also a new solution. As in the 
charged Sultana-Dyer case, the natural source of the electromagnetic field is a central 
electric charge with an additional global electric current. The global structure is briefly 
studied and it is verified that the corresponding spacetime may represent a charged black 
hole in a cosmological background. All the solutions present initial singularities as found 
in the McVittie metric. 
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1. Introduction 


1.1. Preliminary remarks 

A major problem in general relativity and cosmology is related to the effects of the 
cosmological expansion on gravitating local systems, such as the solar system, accompanying 
the expansion of the universe. This question has been the subject of studies since 1933, with 
McVittie pQ. Several models have emerged since then, and what is clear so far is that the 
way the local dynamics is affected strongly depends upon the choice of the specific metric 
to represent the system. Small scale gravitating systems, small compared to the Hubble 
radius, participate in the expansion, but the effect of the expansion being so small that 
can be neglected. In such cases, the global dynamics is well described by the Friedman- 
Lamaitre-Robertson-Walker metric (FLRW), while the local dynamics is independent of the 
cosmological expansion and described by an independent metric. For large scale structures 
compared to the Hubble radius, the cosmic expansion becomes significant and cannot be 
neglected (see, e.g., [2]). In this situation, a fully dynamical metric including the effects 
of the cosmic expansion must be used. In the case of a local electromagnetic system, the 
analysis of [3] shows that a strongly connected system cannot be disturbed by the cosmic 
expansion. Based on this, we can say that in an expanding FLRW universe only weakly 
bounded systems take part in the cosmic expansion [2]. However, it is still not clear what 
the effects of the cosmological expansion are on other, maybe more subtle, features of local 
gravitational systems such as the existence of trapping horizons in a black hole spacetime. 

On the other hand, the description of black hole dynamics is not well settled; some 
points have been set aside as a first approximation to simplify the analysis. Although the 
black hole thermodynamics relies on the discovery of Hawking radiation urn, to avoid further 
difficulties the backreaction has been ignored in many studies about black hole dynamics in 
cosmic backgrounds. More complication arises when the black hole is immersed in a cosmic 
content other than a cosmological constant, as dark matter and dark energy dynamic fields 
(see, e.g., [6] for an analysis of Hawking radiation from black holes asymptotic to the Einstein- 
de Sitter Universe). In this way, if the observations indicate w < —1, the possibility of the Big 
Rip provides another motivation for the study of local dynamics of black holes, because the 
event horizon can follow the expansion and disappear, exposing the central singularity [7j. In 
general, the cosmic content tends to be incorporated into the black hole [3], and neglecting 
the accretion when modeling the black hole evolution results in unrealistic models. 

Recently observational evidence has accumulated, mainly by the observations of SNe la, 
that sets the equation of state parameter asw = — 1.02+Q^g [5], showing that the universe goes 
through an accelerated expansion phase. This fact brought even more interest on the studies 
on the effects of the cosmic expansion on local gravitating systems. Therefore, it would be 
very useful to have exact solutions of the Einstein equations which could describe objects with 
strong gravitational field embedded in an accelerated expanding universe. Further study on 
the already known exact solutions to acquire a full understating of the geometric and physical 
properties of such solutions is also an important task. This is one of the aims of the present 
work. 

1.2. McVittie and Vaidya type solutions 

The first known solution of Einstein equations that has been written as an attempt to describe 
a spherically symmetric pointlike particle in an expanding universe is the McVittie metric [T|. 
In fact, such a metric may be used to investigate the effects of cosmic expansion on local 
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= a\t) (l + 2 ^y) (*'+^V(l-^) ( 1 + *5i)) dt2 ' (1> 


where t is the cosmological time, r is a radial comoving coordinate and a(t) is to be interpreted 
as the cosmological scale factor. It was first though that metric |l]) could describe the 
Schwarzschild black hole geometry in a cosmological background. In this scenario the mass 
to of the point mass (or the central black hole) is constant, there is no flow of energy towards 
the central object, or escaping from it. This non-accretion condition says also that the size 
of the central object does not change during the expansion. However, it has been argued 
(see, e.g., mmm, see also mm) that the McVittie metric cannot describe a black hole 
evolving in an FLRW universe due to the presence of a spacelike singularity at r = m/2a(f), 
which corresponds to an infinite pressure of the cosmic fluid surrounding the object. According 
to these authors, the metric describes a singular central spherical object, but the event horizon 
never develops. A timelike limiting surface and an event horizon are formed, but both shrink 
due to expansion. On the other hand, it is known that in the Einstein-de Sitter universe, in 
which a(t) ~ e Hot , with constant Hg, the McVittie metric reduces to the Schwarzschild-de 
Sitter metric, and the central object described by this metric is a black hole. In fact, more 
recently it has been shown that if the cosmological expansion is such that for late times the 
Hubble parameter H(t) = d(t)/a(t) (with the over dot (') standing for the total derivative) 
reaches a constant value, in which the spacetime asymptotes to the ACDM cosmology, an 
event horizon is formed and black hole and white hole regions are present in McVittie solution 

da Esina esi . 


Straightforwardly from the McVittie metric, one can build new black hole type solutions 
in cosmological backgrounds by replacing the mass parameter to by a function of the cosmic 
time t, to —^ m(f). For instance, by taking m = M a(t ), with constant M, and defining a new 

— j in metric 0 , we find 


radial coordinate R by R = r 1 + 


ds 2 


a 2 {t.) 





( 2 ) 


which is a particular case of the Thakurta metric m- Of course, other choices of m(t) 
yield different geometries that may be interesting as candidates for representing black holes 
immersed in cosmological spacetimes (see below). Since the Einstein equations are not 
invariant under the change to — > m(t) one sees that such a procedure changes the sources. In 
the McVittie solution the source is a perfect fluid with energy density p given by 8np = 3 H 2 (t), 
and pressure p, by 87175 = —3 H 2 (t) — 2 H(t) [1 + m/2ra(t)\ / [1 — m/2ra(t)\. In turn, for the 
Thakurta metric ([2| the simplest source is a fluid with density 8irp = H 2 (t)/ (1 — 2 M/R), 
pressure 8np = — (^3H 2 {t)+ 2H(t)j / (1 — 2 M/R), and with an energy flux Q given by 

47tQ = M H(t)/ R 2 ( 1 — 2 M/R) 3 ^ 2 in the radial direction. Both spacetimes have singularities: 
McVittie at r = m/2a(t ) and Thakurta at R = 2 M. 

With the aim of describing black holes in expanding universes, other metrics have been 
written down. Vaidya i2D: proposed a metric to describe a Kerr black hole immersed in a 
cosmological background, which, in the case with no rotation, can be written in the form 


ds 2 = a 2 (t ) [dr 2 + r 2 dfl 2 ] + 


2 TO 


dt 

, a (t) 


dr 


-dt 2 


r 


(3) 
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where to is a constant. Vaidya argued that the metric represents a black hole whose event 
horizon is located at ra{t) = 2m. However, this seems a more subtle question than initially 
believed, since the surface r = 2m/a(t) does not look like an apparent horizon (see , e.g., 

my 

A similar (Vaidya type) metric as a possible representation of a Kerr black hole embedded 
in the Einstein-de Sitter universe was examined by Thakurta [T9]. As it was shown in [211 , in 
the Vaidya and also in Thakurta solutions some of the energy conditions are violated. Even 
though these possible drawbacks, more studies on these metrics are necessary in order to 
understand the complete geodesic structure of the corresponding spacetimes. 

Another class of solutions found by Sultana and Dyer (SD) [22] is also of interest for the 
present work. The SD metric, 

ds 2 = a 2 (t) [dr 2 + r 2 dfl 2 ] + [dt + a(t)dr] 2 — dt 2 , (4) 

with to being a constant, was proposed to represent a Schwarzschild black hole in an Einstein- 
de Sitter universe. This solution is similar to the McVittie spacetime, but with the no¬ 
accretion condition being released, i.e., in the SD solution there is matter (energy) accretion 
onto (from) the central object. Another difference is that the McVittie solution can be sourced 
by a single perfect fluid, while SD spacetime needs a mixture of two non-interacting perfect 
fluids, a massive dust and a null (lightlike) dust fluid. In SD solution, both in radiation and 
matter dominated era, some of the energy conditions are violated. 

Notice that the SD metric Q can be obtained from the Vaidya metric ([3| by replacing 
the mass parameter to by ma 2 {t). In other words, the Vaidya and the SD metric differ from 
each other only by the rate of mass accretion onto the central body. 

Finally we mention here the attempts made more recently as, for instance, the model 
proposed in 1231 in which the mass of the central object changes due to accretion of a phantom 
fluid, during the cosmic expansion, considering backreaction. This metric differs from the 
solution of [H]; however the results can not be compared effectively since the later does not 
consider backreaction. 

1.3. Charged black holes in expanding universes 

The charged black hole generalization of the McVittie metric was found by Vaidya and Shah 
[231 ; 25]. Later, a generalization of the Vaidya @2] solution including electric charge in the 
source was written |26j . More recently, a very interesting charged solution was put forward 
by Kastor and Traschen pzj. Starting with the Majumdar-Papapetrou multi-black holes 
static solution of [25], the authors found a solution of the Einstein-Maxwell equations which 
describes the dynamics of a system of charged multi-black holes in a cosmological background. 
Among other interesting aspects, the solution furnishes an analytical description of the 
dynamics of coalescing black holes. In the case of a single charged black hole, the solution 
corresponds to the extremal Reissner-Nordstrdm black hole in a spatially flat Einstein-de 
Sitter universe. The global structure of the Kastor-Traschen spacetime was explored in [29] . 
The construction of dynamical multi-centered black hole type solutions was generalized to an 
arbitrary FLRW universe in a recent interesting work (see [30] 3. 

Following a different path, and long before the works on multi-centered charged black 
holes, in [24] ;25j it was put forward a charged version of the McVittie, which is interpreted as 
a charged Reissner-Nordstrom (RN) metric embedded in a FLRW cosmological background. 
The influence of the cosmological evolution on the size of the RN black hole, as well as 
the motion of a test particle close to central charged body, were analyzed in m ■ Several 
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properties of such a spacetime were investigated, but the global structure was not explored 
yet. 

In ST] a set of metrics representing charged black hole type solutions in expanding 
backgrounds was displayed. The authors were interested mostly in verifying the energy 
conditions of the resulting spacetimes. They analyzed different charged metrics. The first 
one was a charged version of the McVittie metric, corresponding to a particular case of the 
solution presented in [25] IE]. A charged Vaidya metric was also considered. Other metrics 
studied there which interest us here were the charged Thakurta m and the charged SD [22] 
metrics. A new solution, similar to the McVittie spacetime, was also considered in [21] ■ The 
energy conditions were studied in all of the solutions and it was shown that some of them 
may be of physical interest. 

Some of the solutions considered in [21] may be obtained from the corresponding charged 
static solutions through conformal transformations, from what some of the properties of the 
resulting spacetime could be inferred. However, other cases such as the Vaidya type metrics 
are not obtained in such a way, and the sources of the electromagnetic held, as well as the 
matter content of the background, must to be analyzed with care. Our aim in this paper is 
firstly investigate such electromagnetic sources in detail, and explore the main properties of 
the corresponding energy-momentum tensor. Moreover, we also briefly investigate the global 
structure of the associated spacetimes. 

1.4. This paper 

The main aim of the present work is to investigate the electromagnetic sources of some 
charged black hole type metrics in expanding universes. Inspired by the charged versions 
of the McVittie [T] and Vaidya [2Dj type metrics previously considered in the literature 
[21ES] El HI], we write a sufficiently general metric with spherical symmetry and show that, 
under certain circumstances, there can be a global electric current in the radial direction. 
Then, some particular cases are studied. In particular, the charged Vaidya, Sultana-Dyer, 
and Thakurta metrics reported in [21] are considered. The mass parameters of these solutions, 
as well as the charge parameters, are taken as a function of the cosmological time, and new 
charged solutions are then built. A brief analysis of each new metric is given. 

The paper is structured as follows. In Section [2] we write a sufficiently general metric 
and give the general form of the Faraday-Maxwell tensor, and of the stress-energy (energy- 
momentum) tensor of the electromagnetic field related to such a metric. The definition of 
apparent horizon is also given in that section. Then, the charged McVittie metric is studied 
in section [5] where attention is paid to the electromagnetic aspects of the solution. Some 
aspects of the global structure of the corresponding spacetime are also considered. The 
studies reported in section [3] are not new, but the results are presented there for comparison 
to the other cases. Section^] is dedicated to study the general properties of a charged metric 
written based on the original work by Vaidya [20] , Here we replace the mass parameter 
of the original metric by an appropriate function of the cosmic time. Again the main 
interest is on the electromagnetic quantities of the solution, and a brief analysis of the global 
structure of the resulting spacetime is given. The results about apparent horizons presented 
in this section are new results. A charged version of a metric originally given by Sultana 
and Dyer [22] is analyzed in section [5j It is shown that the electromagnetic stress-energy 
tensor corresponds to a situation where there is a global electric current across the spacetime. 
Again, the mass and charge parameters of the original solutions are replaced by particular 
functions of the cosmic time and some aspects of the global structure, including the apparent 
horizons, of this new charged SD spacetime are also shown. Similarly, a charged version of the 
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nonrotating Thakurta metric p2j is investigated in Section[6] New mass and charge functions 
are considered, and some aspects of the global structure of the corresponding spacetime are 
also given. In Section [7] we conclude. 

2. A sufficiently general metric, the electromagnetic field and the stress-energy 
tensor 

2.1. The metric, Maxwell equations and solutions 

A sufficiently general metric for all cases we are interested here is of the form 

ds 2 = — fo(r, t) dt 2 + a 2 (t) fi(r, t) dr 2 + 2 a(t) f 2 (r, t)dtdr + a 2 (t ) r 2 f^(r, t) dfl 2 , (5) 

where t is a timelike coordinate, r is a spherical coordinate, a(t) is the expansion factor, dTi 2 
is the metric on the unit sphere, and functions fo(r,t), /i(r, t), f 2 (r,t) and fs(r,t) depend on 
the indicated coordinates. To simplify notation we write gi 2 (r, t ) = /o(r, t ) /i(r, t ) + / 2 2 (r, t). 

Let us stress that we are going to analyze a few particular cases of electrically charged 
spherically symmetric black hole type solutions in expanding universes, for which the functions 
/o, /i, /2 and fs are given. For instance, in the case of the charged Vaidya type metrics 
considered in section |4j one has the simple form fs(r,t) = r 2 . The same holds in the case of 
the SD and Thakurta type metrics considered below. For these three particular metrics, it 
results gi(r, t) = 1, and the radial coordinate r plays the role of a comoving radial coordinate, 
analogously to the comoving radial coordinate in FLRW cosmologies. With this in mind, 
we present in this section an analysis of the electromagnetic fields and sources for the given 
metrics by using the form given in equation (|5j) to represent anyone of these cases. 

The general form of the Faraday-Maxwell tensor F pv and the electromagnetic stress- 
energy tensor for the metric ([5]) are relevant to the present work. As a consequence of 
the spherical symmetry, the only nonzero components of are F tr = —F rt = E(r,t). 
The contribution of the electromagnetic field to the energy-momentum tensor, i.e., the 
electromagnetic stress-energy tensor E v = F pa F va /(An) + 6/F pa F pa / (16n) then is 

El = E r r = -El =-E% = - a ^ g\ 2 (r, t)E 2 (r , t), (6) 

with all the other components of E” being identically zero. The Maxwell equations give the 
two equations 

d ^g r ' ^ = Aira 3 (t)r 2 f 3 (r,t)g 1 (r,t) J*(r,f), (7) 

^ = -Ana 3 (t)r 2 f 3 (r,t.) gi(r,t) J r (r,t), (8) 

where we defined 

Q(r,t) = a 3 (t) r 2 f 3 (r,t) gi(r,t) E(r,t), (9) 

with J 4 (r, t) and J r (r,t) being the only nonzero components of the electromagnetic current- 
density. The other components of the Maxwell equations are identically zero. 

Equation (Jt|) gives that the quantity Q(r,t) = Ana 3 (t ) r 2 fz(r,t ) g\(r,t) T 4 (r,t) dr may 
be interpreted as total electric charge inside a sphere of radius r at time t. This means that, 
using equation Q, we may interpret dQ(r,t)/dt as the electric current across the spherical 
surface of radius r at time t. 
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We are mostly interested in cosmological solutions with no net electric charge distribution 
across the whole spacetime. However, the location of the charged object cannot be determined 
without specifying the explicit form of the metric functions in equation 0 , and, moreover, 
since there must be some electric charge somewhere in the spacetime, the nonstationary 
character of the metric implies there exist some charge flux, i.e., a current-density which, 
of course, depends on the choice of the coordinates. If one takes J*(r, t) = 0, one has 
from equation 0 Q(r,t) = Qft ), a function of time alone. The radial component of the 
current-density J r (r,t) may be nonzero. Therefore, there are at least two different solutions 
to equations Q and (|8]) which are interesting for the present analysis. One of which holds 
in a region of the spacetime where there are no electromagnetic sources, J 4 (r, t) = 0 and 
J r (r,t ) = 0, while the second is when electric current is allowed, i.e., for a special choice of 
coordinates, one may choose J 4 (r, t) = 0 but with J r (r,t ) ^ 0. In the next sections, when 
studying charged versions of McVittie, Vaidya, SD and Thakurta metrics, we return to this 
subject. 

We can determine the total electric charge inside a spherical surface at a fixed radial 
coordinate r for the spacetimes represented by metric |5]). Let S represent the spherical 
surface (t, r) = constants. The electric charge may be defined as 

Q{r,t)= f F^n v dS, (10) 

Js 

where C is the unit normal to surfaces of constant t, n v is the spacelike unit vector orthogonal 
to S pointing outwards, and dS is the area element on S. Vectors and n M also satisfy the 
relations Ct^ = — 1, t^n^ = 0, and nTn^ = 1. In the case of metric ([5]), we may choose 
^ = 8tVfi(r>t)/gi(r, t ) - Sff 2 (r, t)/[a(t)gifr, t)y //i(r, t)} and nT = 5p / yjfi(r,t), and the 
resulting expression coincides with equation ^. Such a definition of electric charge is applied 
to some particular cases below. 


2.2. The scale factor a(t) 

For the charged McVittie, Vaidya, SD and Thakurta type metrics which are considered in 
the present work, the function aft) in equation ([5]) corresponds to the cosmological expansion 
factor. In general, this function is defined by the matter content of the universe, after solving 
the full Einstein-Maxwell system of equations. However, regarding the kind of study being 
performed here we satisfy ourselves by taking two particularly simple cases, defined below. 
The first choice is, for comparison, the expansion factor used in nn nu 


a(t) 



CH) 


with k being a constant to be chosen appropriately. This choice has an initial power-law 
expansion with aft) ~ f 2 / 3 and a final de Sitter accelerated phase. From now on this choice 
is called the case (a). 

The second choice is a power-law expansion, 


a(t) = i“, (12) 

with constant a, which corresponds to a universe filled with a perfect fluid whose equation of 
state is of the form p = lo p, with uj = (2 — 3a)/2. Our choice shall be a = 2/3, and so w = 0, 
a cold dark matter dominated model. This is called the case (b). 
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2.3. Singularities and horizons 

Besides the study on the electromagnetic sources, we investigate the curvature singularities 
and horizons of the chosen metrics. 

The interest on curvature singularities in the cosmological scenario is because they are 
physical, i.e., true spacetimes singularities, meaning that physical quantities such as the energy 
density or pressure of the cosmological fluid diverge at the singularity. It is then necessary 
to calculate the usual curvature invariants which are built from the Riemann tensor R pvpa . 
As argued in BE for a metric of the form given by equation ([5]), the Ricci (TZ) and the 
Kretschmann (1C) scalars are necessary to identify the spacetime singularities. These are 
calculated for each metric investigated in the present work (see below). 

In order to analyze the global properties of a spacetime whose metric has the form of 
equation ©, it is useful to perform a coordinate transformation which brings in a new radial 
coordinate R defined by 

R = a(t) r \J fo(r, t). (13) 

With this, metric ([5]) assumes the form 


ds 2 


h(R , t) 


C 2 (R,t)' 
A(R,t ) _ 


dt 2 + A(R, t ) dR 2 


2 C(R,t)dRdt + R 2 dfl 2 , 


(14) 


where the functions A(R,t ), h(R,t ) and C(R,t) are given in terms of the original functions 
fo(r,t), fi(r,t), / 2 (r, t) and fs(r,t), but we do not write such relations here. 

At this point it is interesting to find the equations for the expansions of the outgoing 
(+) and ingoing (-) congruences of null geodesics, 6± = V ;i /c±, where k± are tangent vectors 
to the radial null geodesics. Using this definition and the geodesic equations for the metric 
dl4l) , it follows 


i± = 


R 


C(R,t) 

A(r,t) 


± 


I h(R,t) 
A(R,t ) 


(15) 


dt 


where — ) = kh, p being an affine parameter along the geodesic curve. For well behaved 

\dpj ± 

k l ± , one then sees that the zeros of 6± are given by the zeros of the function between brackets 


in equation (15). This means that surfaces for which the product 6 + 6_ = 0, if they exist, are 


located at regions of the spacetime satisfying the equation 


Fh(R, t) = 


h(R,t ) / C 2 (R,t) 


A(R,t) \h(R,t) A(R,t) 


- 1 = 0 . 


(16) 


According to Nolan BE which follows [32j . for a spherically symmetric metric, the zeros of 
Fh(R, t), which are in fact the zeros of the function y (i?, t) = (V M i?) V M R, defines a trapping 
boundary. For the kind of metrics we are interested here, these are also trapping horizons. On 
the other hand, other authors (see, e.g., pans Ei) call such surfaces generically as apparent 
horizons. We follow the nomenclature of mm- 

We can also verify that the zeros of the expansions 6± coincide with stationary points 


(regions) of the radial null geodesics of the metric (14). In fact, the radial null geodesics are 
given by solutions of the equation 


dR\ 

dt )± 


1 


A(R,t ) 


(c(R,t) ± y/h(R,t)A(R,t)) , 


(17) 
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where the plus (minus) sign indicates outgoing (ingoing) geodesics. We see also that the 
functions ( dR/dt) ± change sign at the zeros of the function T#(I?,f), meaning that the zeros 
may be apparent horizons separating the spacetime into trapped and non-trapped (regular) 
regions. 

It is worth mentioning that, for the metric of equation (14), one has Fn(R,t) = 
h(R,t)g RR , with g RR being the component of the inverse metric. Indeed, this procedure of 
looking for the zeros of the function g RR to find the possible apparent horizons of black hole 
type metrics in expanding spacetimes was followed by many authors (see, e.g., mum nn ng 
for McVittie spacetime, for Kastor-Traschen solution E3, and m for a Thakurta type 
metric). 


2-4- About the next sections 

In the following we investigate a few particular charged metrics written to describe charged 
black holes in expanding universes. The main aim is to understand the physical nature of 
the electromagnetic sources of the corresponding spacetimes. Another interest is to look for 
apparent horizons in each case, the ultimate goal is to determine if black hole horizons are 
really present, besides finding the spacetime singularities in order to identify the physical 
nature of the central object in each case. In order to do that, the complete geodesic structure 
of the given spacetime must be investigated. However, this study is lengthy and, moreover, 
all the metrics we consider here deserve a separate detailed analysis because of the interest on 
the subject in recent literature. Hence, to avoid a very long paper, in this work we investigate 
only the form of the apparent horizons and the curvature singularities for each one of the 
studied metrics. 


3. Charged McVittie metric 

The generalized version of McVittie metric including the electric charge of the central body 
was put forward in [25] (see also [31]). Here we review the electromagnetic sources of such a 
metric, and present a short analysis of curvature singularities and apparent horizons. 

3.1. The solution 

In the case of an asymptotically spatially flat FLRW spacetime, the charged McVittie 
metric may be obtained from the general form © by setting f 0 (r,t) = f 2 {r,t)/g 2 (r,t), 
fi(r,t) = fs{r,t) = g 2 (r,t), and f 2 (r,t) = 0, which gives gi(r,t) = f(r,t), where the functions 
f(r,t ) and g{r,t) are defined respectively by 


= 

1 m + q 

4 a 2 (t)r 2 4 a 2 (t)r 2 

(18) 

9{r,t) = 

( m \ 2 q 2 

\ 2 a(t)r) 4a 2 (t) r 2 ’ 

(19) 


with m and q being constants related respectively to the mass and charge of the central body. 
The resulting metric is 

ds 2 = — ^ ^ ^ dt 2 + a 2 (t)g 2 (r, t) (dr 2 + r 2 dtt 2 ) . ( 20 ) 

g 2 (r,t) 

The function a(t) is interpreted as the expansion factor with t being a cosmological time. 
The asymptotic metric (r — > oo) results ds 2 = — dt 2 + a(t ) 2 ( dx 2 + x 2 dfl 2 ), which is the 
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FLRW metric in the case of flat three space. The above metric suffers from the same illness 
as the uncharged McVittie metric, but at least when the Hubble parameter H(t) = ait)/a 
asymptotes a constant at t —> oo, such a solution has an event horizon and represents a black 
hole in an expanding universe (see, e.g., nmng). 


Notice that a further generalization of the charged McVittie metric (201 is obtained 


simply by replacing the constant parameters m and q by function of the cosmological time, 
m(t) and q(t). As discussed in the literature (see e.g. [TUllTTllT3] h this implies in the violation 
of the non-accretion hypothesis of McVittie, and there must be some kind of energy flux in 
the radial direction throughout the spacetime. In respect to the electric charge, it means 
a non-constant electric charge of the source, and the presence of a radial electric current 
throughout the spacetime (see also the next section). 


3.2. Faraday and stress-energy tensors, conserved electric charge and current-density 


The electromagnetic source and the stress-energy tensor for the metric (20) were analyzed 


in mm and then we do not need to reproduce the results here. However, further comments 
with respect to the Maxwell equations are in order. According to our analysis in section[2j the 
nonzero components of the Faraday-Maxwell strength tensor are F tr = —F rt = E(r,t), and 

for the charged McVittie metric one finds (see also [2HI EH) E(r,t) = —sv—, ■ 

a 6 (t) r 2 f(r, t) g-{r , t) 

The solution presented in [25] assumes h(t) = 1, leading to Q(r,t) = qo = constant, and go is 
a conserved quantity identified with the total electric charge of the source, so that qo = q, q 
being the charge parameter of the metric. This result means that comoving observers in the 


charged McVittie spacetime (20) attribute a constant electric charge to the central object. 


Furthermore, with such a choice, the Einstein-Maxwell equations in the presence of a cosmic 
perfect fluid are satisfied resulting in Friedmann type equations for the energy-density and 
pressure. No additional terms in the density or pressure arise because of the presence of the 
electric charge. 

On the other hand, if one chooses h(t) ^ 1, the Maxwell equations are satisfied with 

qo h(t) 


the inclusion of a radial current-density given by J r (r, t) = — 


and the 


4Trr 2 a 3 (t) g 2 {r, t) ’ 

charge of the central object is not constant. In addition, in such a case there would be 
additional contributions to the energy density and pressure of the cosmic fluid coming 
from the electromagnetic field. For instance, the energy density has a term of the form 
q 2 (t)/i A — qQh 2 (t)/r A , which, for h(t) ^ 1 and constant q implies in violation of the energy 
conditions. This solution has not been investigated in the literature, since it is regarded as 
unphysical. It is worth mentioning that if we replace the parameter q by a function of time 
and take q(t) = qoh(t), the Maxwell equations are satisfied and the spurious term in the 
energy (and pressure) of the perfect fluid vanishes. However, in order to satisfy the Einstein 
equations additional matter components must be added to the background. A detailed study 
of such a case is not the purpose of this work. 


3.3. Singularities and horizons 

As already mentioned, we do not investigate the complete geodesic structure of the charged 
McVittie metric here, since a complete analysis is lengthy and out of the scope of this paper. 
However, the results on the singularities and horizons are of interest. We consider here just 
the case of constant mass and charge parameters m and q , and show the singularities and 


apparent horizons of metric (20) for comparison to the other solutions investigated in the next 
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sections. More details about the apparent horizons and singularities in the case of constant 
parameters m and q can be found in [ 551 . 

To say something about curvature singularities of the charged McVittie metric (201, we 
calculate the Ricci and the Kretschmann scalars (see Appendix AI. To simplify expressions 


it is convenient to use a new radial coordinate R defined by R = a{t)r g(r, t) (see fBT]). The 
Ricci scalar is then 1Z = 12 H 2 (t) + 6H(t)/y/l — 2 M/R + Q 2 /R 2 , where we have written 
M = m and Q = q. The Kretschmann scalar is written in |Appcndix A| We immediately 
see possible singularities when H(t) diverges, and at the region of the spacetime where 
h(R) = 1 - 2 M/R + Q 2 /R 2 = 0, i.e., for R = M± yjM 2 - Q 2 , or, in terms of the original 
coordinate r, at r a{t) = ±s/M 2 — Q 2 / 2. In the overcharged case (M 2 < Q 2 ), the function 
h(R) has no real roots and the singularity is at R = 0 (at r a(t) = — (M ±Q)/2). We consider 
here just the undercharged case ( M 2 > Q 2 ). 

The apparent horizons of the McVittie metric are given by the roots of the equation (see 
equation ©) 

2 M Q 2 ' 


F M {R,t) = ~H 2 (t) R 2 + 




( 21 ) 


which is a fourth-order polynomial equation for R, given the solutions in terms of the time 
t. Note also that the zeros of F M (R : t) coincide with the zeros of the radial coordinate 

velocity dR/dt of ingoing lightlike geodesics, i.e., the zeros of ^ = H(t) R\J 1 — ^ ^ — 

(l “ + w) > for which 1 - 2 M/R + Q 2 /R 2 ± 0. 



Figure 1. The evolution of the apparent horizons as a function of the time t in the 
undercharged case of the charged McVittie metric. We used M = 1.5, Q = 1.0. The case 
(a) is for the scale factor of equation 0 with k = 0.1, and the case (b) for a(t) given by 
equation ( | 12 [ ) with a = 2/3. For sufficiently large times there are three real positive roots 
of the equation FM(R,t) = 0, with Fm(R, t) given by equation ( |21| ). There is an initial 
singularity at R = R+ = M + yjM 2 — Q 2 , t = finite, Ri(t) is always smaller than R .|_, 
l? 3 (t) is a cosmological type horizon while R2(t) is a black hole type apparent horizon. 


For the undercharged McVittie metric ( M 2 > Q 2 ), and for the chosen form and 
parameters of the scale factor, for which k x M is sufficiently small (case (a)), and a = 2/3 
(case (b)), function Fm{R , t) has just one real root at small times, while it has three real roots, 
Ri(t), R 2 (t) and R${t) for sufficiently large times (see figure [l]). The root R\(t) is always 
inside the singularity radius R + = M + \JM 2 — Q 2 and so it is of no interest, as long as we 
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are interested in the R > i?+ region. In view of the similarity between functions R 2 (t ) and 
I ?3 (t) shown in figure [I] and the corresponding functions of the uncharged McVittie case, we 
may infer that the geodesic structure of the undercharged McVittie is similar to the original 
McVittie spacetime, at least for R > R + (see e.g. [TUI Ell HU HSJ EH EH] ) • If the cosmological 
expansion is such that lim H(t) = Hq =constant (case (a)), there are two horizons in the 

t—¥ OO 

future of the initial singularity. They correspond to the spherical surfaces R c and Rh which 
are the asymptotic values of the functions R 3 (t) and Riit) at t —> oo, respectively. According 
to the study of l..(j . the black hole horizon TL is defined by (i? 2 , t = oo), while R c defined 
by (R 3 , t — > oo) gives a de Sitter horizon (case (a)), or the Hubble horizon (case (b)). Our 
Rh here corresponds to the asymptotic limit of R- of nu- The apparent horizon Ri(t) is 
inside the singularity, but it has an interesting interpretation, since it asymptotes a charged 
black hole inner (Cauchy) horizon at t — > oo, i.e., the region (J?i, t = oo) is similar to the 
Reissner-Nordstrom-de Sitter inner horizon; it is a Cauchy horizon. In case (b), in which 
lim H(t) = 0, the conformal structure of the boundary of the charged McVittie spacetime is 

t—> OO 

similar to that of a Reissner-Nordstrom black hole, apart from the naked initial singularities. 
In this case, the asymptotic form of the metric in (R,t) coordinates is time independent. In 
both cases (a) and (b), the asymptotic metric, at t —> oo, can be put into stationary forms, 
so that both admit Killing vectors associated to time translations. It is interesting to explore 
the geodesic structure of the charged McVittie spacetime in detail, including the extremal 
and overcharged cases. However we are not going to report the results of such a study here. 


4. Charged Vaidya type metrics 

A charged version of the Vaidya metric [20| was reported in ED- Here we review the 
electromagnetic sources of a modified version of that metric, and present a brief analysis 
of curvature singularities and apparent horizons. The modifications of the metric is done by 
considering also time-varying mass and charge parameters when compared to the solution 
given in [2T1 . 


4-1- The solution 

According to our notation of section [2j the charged Vaidya type solution given in m is 
obtained by choosing f 0 (r,t) = 1 - g{r)/a 2 {t), fi{r,t) = 1 + g(r)/a 2 (t), f 2 (r,t) = g{r)/a 2 (t), 
and f 3 (r,t ) = 1, with g[r) = 2 m/r — g 2 /r 2 , which gives gi(r,t) = 1 (see equation ([5])). Here 
we consider a generalized version of the metric given in [21] by taking the mass and charge 
parameters as a function of time, i.e., to = m(t) and q = q{t ), and so we write 


ds 2 = a 2 {t) \dr 2 + r 2 dQ 2 ] + g(r, t) 


dt 


where 


g(r,t) = 


2m(t) 


a(t) 
q 2 (t) 


+ dr 


-dt 2 


( 22 ) 


(23) 


As promptly seen from (22), for a(t) = 1 and constant to and q it results the Reissner- 
Nordstrom black hole spacetime, and far from the central body, r —> oo, it results the Einstein- 
de Sitter cosmological spacetime. 

The nonzero components of the Einstein tensor corresponding to metric (22) are 
the diagonal ones G t t , G r ri G e g = G v besides the ones that mix the radial and time 
coordinates, namely G t r and G r t ■ We do not write the full expressions of these components 
here, but we observe that there are changes in the sources in comparison to the results given in 
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[2l ! because of the time dependence of m{t) and q(t). To be explicit, we take, for instance, the 
new terms that arise in the component G t t , which are —2 H((t) ( m(t)a{t)r — q{t)q(t)) /a 2 (t)r 2 . 
As in the original Vaidya metric [20] . these terms are singular at r = 0. Their contributions to 
the energy density depend on the explicit form of the functions rn(t) and q(t ), and in the case 
m(t) = M/a(t), q(t) = Q/a(t), considered in the next subsection, the problem of negative 
energy densi ties (at early t imes t and large ?’) of the original Vaidya solution persists (see, e.g., 
[2T1 . see also Appendix A). We need to point out explicitly the linear terms on q 2 (t) and that 
do not contain derivatives of a(t), because they are directly related to the electromagnetic 
stress-energy tensor, what is one of our main points of interest. The contribution of these 
terms to the Einstein tensor are the following 


G t _ /-rr _ /^i( 

t — Lr r — —Lx 


q 2 (t) 


1 _ ___ 

v ~ a 4 (f)r 4 ’ 


(24) 


From the analysis of section [2] it is expected that terms of the form q 2 (t)/a‘ i (t) r 4 would be 
exactly canceled by the corresponding components of the electromagnetic stress-energy tensor 


E”. This fact is confirmed in the next subsection. 


4-2. Faraday and stress-energy tensors, conserved electric charge and current-density 

Following the analysis of section [2j one gets the Faraday-Maxwell tensor field as E(r , t ) = 
qoh(t) /a 3 it)r 2 , qo being a constant, which yields Q(r,t ) = qoh[t) (see equation 0) , and the 
electromagnetic stress-energy tensor then is 


jpt _ zpr _ rpv 

Fj t — 


= _E>!P = ~go fe2 (^) 

v 87 t a 4 (f)r 4 


(25) 


with the other components being identically zero. 

It is seen that by choosing h(t) = 1, i.e., E(r,t) = qo/a 3 {t)r 2 , with constant q 0 , the 
Maxwell equations are satisfied with zero current-density, J 4 (r, t) = 0, which implies that 
q(t) = qo is the total conserved electric charge of the source, as expected. In fact, from Maxwell 
equations 0 and 0 , we see that the electric four-current is not well defined just in the region 
of the spacetime for which a 3 (t) r 2 = 0. Everywhere else in the spacetime the electromagnetic 
sources vanish, and then it is natural to associate the loci where a 3 {t)r 2 = 0 as the actual 
location of the source. Furthermore, comparing equations (25) with equations (24), it is seen 


that the components of the electromagnetic stress-energy tensor cancel the corresponding 
terms of the Einstein tensor through Einstein equations ( G = 87 t E^). 

Contrarily, if one takes h{t) 7^ 1, so that E(r,t) = qoh{t)/a 3 (t)r 2 , the radial electric 

current-density is nonzero, J r (r,t ) = — . ^ . , and the total charge of the source is 

47 t r-a 6 (t) 

not a constant. In this case, there may arise contributions to the energy-momentum tensor 
from the electromagnetic field that violate the energy conditions. In fact, similarly to the 
uncharged case (as investigated, e.g., in [2D]), the charged Vaidya metric (22) can be sourced 
by a global perfect cosmic fluid in addition to a heat flux, besides the spherical charged 
object. If h(t) 7^ 1, the energy density (and also the pressure) of the cosmic fluid acquires 
a term proportional to q 2 (t) — q^h 2 (t). The time component of the Einstein equations gives 
Gj — 87 tE\ = — ( q 2 (t ) — q 2 h 2 (t)) /r 4 a 4 (t). Th is is similar to what happens in the charged 
McVittie metric we investigated in subsection 3.2 if q 2 (t) 7^ <7gh 2 (f) the energy conditions 
may be violated and the corresponding solutions is disregarded as unphysical. 
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Let us comment here on the curvature singularities and possible apparent horizons of the 
spacetime defined by the charged Vaidya type metric (221. 

In the case we choose constant m and q , the Ricci scalar is given by 7 Z = 
2 H 2 (t) (6 + g(r)/ a 2 (t)) +2 H(t ) (3 + 2 g(r)/a 2 (t)) — AmH(t) /a 3 (t)r 2 , while the Kretschmann 
scalar is given in Appendix A| We can note that both scalars are singular when H(t ) —> oo, 
and at a(t) r —> 0. 



Figure 2. The evolution of the apparent horizons as a function of the time t in the 
undercharged case of the charged Vaidya metric. We used M = 2.0, Q = 1.0. The case 
(a) is for the scale factor of equation O with k = 0.2, and the case (b) for a(t) given by 
equation ( |12| ) with a = 2/3. For sufficiently large times there is only one real positive root 
of the equation Fy(R,t) = 0, with Fy(R,t ) given by equation ( |26| ), which corresponds to 
a cosmological (Hubble like) horizon. There is an initial singularity when o.(l ') =0 and 
also at R = 0, t = finite. 


For the non-constant m(t) and q(t) case there are of course many interesting choices. To 
keep the metric similar to the McVittie one, and to compare with the SD and Thakurta type 
metrics in the following, we take the case m(t) = M/a(t) and q(t) = Q/a(t) with constant M 
and Q. A long but straightforward calculation shows that the curvature scalars are singular 
at the same regions as in the case with constant m and g, where a(t)r = 0. As commented 
in Appendix A the curvature scalars show that the singularities are the same as for constant 


m and q. 

Similarly to the case of the charged McVittie metric, the trapping horizons of a spacetime 
whose metric has the form of equation (22) may be found by analyzing the roots of the 
equation (see equation ©), 


F V (R, t) = H\t ) R 2 + [H{t) R - 1]' 


a 2 (t) 


-1 = 0 , 


(26) 


where now R = ra(t), g(R) = 2 M/R — Q 2 /i? 2 , and we defined M = ma(t) and Q = qa(t). 
Here we consider the case of constant M and Q only. 

The zeros of the function Fy(R , t ) tell us where the expansion rate of radial null geodesics 
changes signs, which is a necessary condition to find apparent horizons. The short analysis 
we give here is by comparison to the McVittie metric studied in mm- The solutions of 
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equation (261 depend explicitly on a(t) and hence we need to define such a function a priori. 
In the present analysis we use the two cases mentioned in section [2. 2 1 


Let us first mention that equation (26) has at least one real positive root for all times 
given by H(t)R = 1. This solution is shown by the curves R^t) in figure [2] it is the 
Hubble radius, it tends asymptotically to a cosmological de Sitter horizon in case (a), with 
R = l/H(oo) = constant, and is exactly 3£/2 in case (b). Moreover, since in an expanding 
cosmological scenario the scale factor becomes arbitrarily large for infinite time, the second 
term in equation (26) vanishes and i? 3 (£) is the only root at late times. For simplicity we 
consider here the undercharged case only, M 2 > Q 2 . The roots Ri(t), R 2 (t), and i? 3 (t), as 
function of the time t, are plotted in figure [2] for the case with M = 2.0 and Q = 1.0, and for 
the scale factor a(t) labeled as cases (a) and (b); see section 2.2 Since the asymptotic form (for 
infinite time) of the roots of the function Fy{R,t) dictates the structure of the boundaries 
of the given spacetime, we see that a black hole horizon analogous to TL = ( R 2 , t = oo) 
of the charged McVittie metric never forms in the charged Vaidya metric. The black hole 
type apparent horizon R 2 (t) is formed after the initial singularity (at the same time as the 
inner apparent horizon i?i(f)), expands faster than R\(t) at early times, but starts to shrink 
after a given time while R\{t) keeps growing slowly with time. Finally, after a finite time 
interval, both apparent horizons coalesce and disappear. The singularity at R = 0 is naked 
for all later times. Assuming, for instance, a matter dominated universe at late times, then 
Hindoo a(t) H(t) = 0 (case (b)) the only real positive root increases indefinitely with time, 
corresponding to the Hubble radius. When the asymptotic form of the scale factor a{t) 
corresponds to an accelerated expanding universe (case (a)), where lim^oo H(t) = Hq= 
constant, the only real positive root tends to a constant, which corresponds to a de Sitter 
cosmological horizon. Even though these partial results have been found, a more detailed 
study is necessary to determine the global structure of the spacetime generated by metric 


( 22 ). 


5. Charged Sultana-Dyer type metrics 
5.1. The solution 

A charged version of the SD metric, which is obtained through a conformal transformation 
of the Reissner-Nordstrom metric written in Eddington-Finkelstein coordinates, was given in 
[21J . For the present analysis we take the following generalized SD type metric 


ds 2 = a 2 {t) [dr 2 + r 2 dfl 2 ] + g{r,t) [dt + a{t) dr} 2 — dt 2 


(27) 


where again we have put g(r,t) = 2m(t)/r — q 2 {t)/r 2 (see equation As noted 

earlier, the factor a(t) is interpreted as the cosmological scale factor, with t being the 
cosmological time. This charged SD metric is obtained from our general metric ([5]) through 
the identifications fo{r,t) = 1 — 2 m{t)/r + q 2 (t)/r 2 , fi(r,t) = 1 + 2 m(t)/r — q 2 (t)/r 2 , 
f 2 {r,t) = g(r,t) = 2 m(t)/r — q 2 (t)/r 2 , and fsir.t) = 1; from what follows gi(r,t) = 1. 


One can see that the SD metric (271 follows from the Vaidya metric by substituting m(t ) and 
q(t) respectively by m(t)a 2 (t) and q(t) a(t) into the metric (22). 

The uncharged version of metric (27) has been considered in the literature as a good 


candidate to represent a Schwarzscliild-like black hole immersed in an expanding cosmological 
background (see e.g. EH S3]). On the other hand, this charged version (|27| presents an 


additional issue that deserves to be considered carefully. For future reference let us note that 
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metric (271 yields an Einstein tensor whose relevant components presents terms like 

q 2 (t) 


G t fir fiu _ fitp 

t — Lr r — —Ur Q — —Ur <0 —- 


2 (t) 


r 4 ' 


(28) 


Compared to the Vaidya metric case, equation (24), the difference is just by a factor of a 2 {t). 


In the following subsection we give the Faraday-Maxwell and the electromagnetic stress-energy 


tensors for the charged SD metric and compare to Einstein tensor (28 1 to analyze the source 
of the electromagnetic field. 


5.2. Faraday and stress-energy tensors, conserved electric charge and current-density 

According to the study of section [2j the Faraday-Maxwell tensor field for the SD metric ([4]) 
is of the form F tr = E{r,t ) = qoh(t) /a 3 {t)r 2 , where qo is a constant and h(t) is an arbitrary 
function of time, and then the electric charge obtained from equation © is Q(r,t) = q 0 h(t). 
The nonzero components of the electromagnetic stress-energy tensor is 


TT't _ j7ir _ rpt 

JCj+ — -C/~ — — JCJf 


= = -go fe2 (*) 

^ 8ira 4 (t)r 4 


(29) 


Comparing this with the Einstein tensor components (equations ( |28| )) we get terms of the 
form G\ — 87 tE\ = ( qo 2 h 2 {t) — q 2 (t) a 2 (t )) / (a 4 (f)r 4 ). The exact cancellation of such terms 
is expected, and then we must have q 2 {t) = q$h 2 (t) / a 2 (t). 

If one takes q(t) = qo = constant, as the solution given in [ 2 Tj . the expected cancellation 
happens only if h(t) = a{t). However, in that case the constant q 0 = q is not the conserved 


electric charge. In fact, using the charge definition (10) we get Q(r,t) = qoa(t). If so, there 
would be a nonzero electric current in the radial direction J r (r,t) = —qoH(t)/{ATra 2 (t)r 2 ), 
implying that the total charge of the central object varies with time. This is the solution 
considered in [ 21 ], but the presence of the radial electric current was not noticed there. 

On the other hand, the natural choice is to take Q(t ) = qo = constant, where the total 
charge of the central body is constant (cf. equation (|l 0 |)), with h(t) = 1 , which corresponds 
to the case of zero four-current density throughout the spacetime, = 0. In such a case, 
one must have q 2 {t) = < 7 §/a 2 (i), giving rise to a new metric of SD type. To see how this 
choice affects the sources let us take the full expression for G\ — 8nEl in the case with 
m(t) = M/a(t) and q(t) = Q/a(t), where M and Q are constant parameters. This gives 

AM Q 2 \ , AMH(t) 


G\ — 8-kE\ = —H 2 (t) 


3 + 


a(t)r a 2 (t)r 2 


+ 


a 2 (t)r 2 


In the original charged SD metric of 


pi ; one has — 87 rE* = —3 H 2 (t) ( 1 + 


2m 


Hence, apart from the time 


4 mH(t) 
r r* / a{t)r 2 

dependence of the mass and charge functions, this component of the Einstein equations has 
no significant differences in comparison with the solution given in |21| . Similar changes occur 
in the other components of the Einstein tensor, but we do not analyze them here. 

We then conclude that the simplest case with q(t) = qo = constant and equal to the 
electric charge of the spherical source is physically not interesting. The Einstein-Maxwell 
equations are satisfied only with the presence of an additional energy-momentum terms that 
may violate the energy conditions. The energy density acquires an extra term of the form 
q 2 (l — a 2 {t )) / (a 4 (f)r 4 ). This term changes sign at time t such that a 2 {t) = 1, which is an 
arbitrary time, since the scale factor a(t) can be normalized in such a way that a 2 (t) = 1 at 
any given time. A similar situation happens in the case of the charged Thakurta spacetime, 
which we investigate in more detail below. 
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Here we comment briefly on the causal structure of the charged SD spacetime given by the 


metric (27). The Ricci and the Kretschmann scalars for this metric are written in Appendix 
p\j first for constant m and q and afterwards for m(t) = M/a(t ) and q(t) = Q/a(t), with 
constant M and Q. As in the case of charged Vaidya solution, both scalars are singular when 
H(t) —i oo and at a(t) r = 0. The singularities are the same when the mass m and the charge 
q depend on time, as it can be seen by calculating the Ricci and the Kretschmann scalars 
for such a case. Again, we consider here the case where m(t) = M/a(t) and q(t) = Q/a(t). 
As it is seen below, this choice implies the metric Q has a geodesic structure similar to the 
charged McVittie metric. 
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Figure 3. The evolution of the apparent horizons as a function of the time t in the 
undercharged case of the charged Sultana-Dyer metric \27\. We used M = 1.5, Q = 1.0. 
The curves (a) are for the scale fact or g iven by equation ( |11[ > with k = 0.2, and the curves 
(b) are for a(t) given by equation with a = 2/3. In this case, there is a singularity 
at R = 0, t = finite. For late times, th ere are three real positive roots of the equation 
Fso(R,t) = 0, with FgE>(R,t) given by ( |30| l. 


Looking for singularities in metric (27) we see that the locus where gtt = g rr = 0 coincides 
with the horizons of the Reissner-Nordstrom metric, r = r±(t) = m(t ) ± \Jm 2 (t) — q 2 (t), 
which are roots of the function f(r,t) = 1 — 2 m(t)/r + q 2 (t)/r 2 . However, these are neither 
spacetime singularities nor apparent horizons. To find the locus of an apparent horizon, we 
proceed as in the case of the preceding section and transform to a new radial coordinate given 
by R = ra{t). Again, apparent horizons can be found by analyzing the roots of the equation 


g RR (R,t) = 0. In the case of SD metric (27), this gives 


FsoiR, t) = H 2 {t)R 2 + [1 - H(t) R} 2 g{R, t) - 1 = 0, 


(30) 


where now g(R,t ) = 2 M/R — Q 2 /R 2 , with M = m(t) a(t ) and Q = q(t) a(t). The resulting 
equation to be solved is a fourth-order polynomial equation whose zeros depend on H(t), Ad 
and Q. For simplicity, we assume here the condition M 2 > Q 2 , with Q / 0 and consider 
the two cases for a(t) as earlier, i.e., a(t) given by equations (11) and For a given 

range of the parameters M and Q , and for the two given choices of a(t), the real roots of 
equation (30), Ri{t) (i = 1, 2, ...), result in the curves drawn in figure [3j At early times, 


there is one real positive root i? 3 (t) which is a cosmological apparent horizon (or the Hubble 
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radius) and then the initial singularity at R = ra(t) = 0 is not hidden by any kind of horizon. 
Indeed, the root 7? 3 (f) is exactly the Hubble time, i.e., the inverse of the Hubble parameter, 
I? 3 (t) = 1 /H{t). At infinitely large times Rz(t) tends to a constant (7? 3 (t) —)• 1/k) and will 
form the de Sitter cosmological horizon in case (a), i.e., (R 3 , t = 00) is a de Sitter horizon, 
or it will be the Hubble surface (I? 3 (f) = 3t/2) for late times in case (b). For sufficiently 
large times, two other solutions Ri(t) and I? 2 (i) become real and positive. The three roots 
exist for all later times and, interestingly, their asymptotic form suggest the formation of 
charged black hole and white hole regions, with a cosmological horizon. For instance, if 
lim^oo H(t) = H 0 =constant, case (a), the three roots tend to fixed values Ri(t) —>• I?i c , 
/?, 2 (£) —I R 2 c, and Rs(t) —>■ R^ c = 1/k. We claim that these values correspond respectively to 
a Cauchy horizon, a charged black hole horizon and a de Sitter (cosmological) horizon. The 
de Sitter horizon is defined by (R 3 , t = 00), the black hole horizon R is given by (R 2 , t = 00), 
and the inner horizon is defined by (I?i, t = 00). In the case of matter dominated universe, 
case (b), Hindoo H(t) = 0, R\(t) and I? 2 (i) also tend to constant values with properties of 
Cauchy and black hole horizons, respectively, while i? 3 (f) tends to 3t/2, the Hubble radius. 
In both cases, the charged SD metric (27) seems to represent a charged black hole in an 
expanding universe. 


6. Charged Thakurta type metrics 


6.1. The solution 

The original Thakurta m solution is conformal to the Kerr metric and was put forward to 
represent a rotating pointlike mass in an expanding universe. For the non-rotating case, the 
Thakurta metric is conformal to the Schwarzschild metric and is also a candidate to describe 
the gravitational field of a point mass in an expanding universe. The resulting metric in this 
case was analyzed in [ML even though the name used there was the SD metric. A charged 
version of the non-rotating Thakurta solution was given in [21] . Such a metric is of the form 

ds 2 = —f(r, t)dt 2 + a 2 (t) [/ _ 1 (r, t)dr 2 + r 2 dfl 2 ] , (31) 

where 


/M) = 1 - 


2 m{t) q 2 {t) 


(32) 


m(t) and q(t) being related respectively to the mass and electric charge of the central body, 
and are fixed parameters in the original solution. Metric (31) may be obtained from our 
general metric © by choosing f 0 (r,t) = fi 1 (r,t) = 1 - 2m{t)/r + q 2 {t)/r 2 , f 2 (r,t ) = 0 , and 
/ 3 (r, t) = 1, so that gi(r, t ) = 1. As mentioned above, if m and q are both constant parameters, 


metric (31) is conformal to the Reissner-Nordstrom metric, as seen by redefining the time 
coordinate by writing dt = a(r) dr. Despite having been briefly analyzed in [21] OB] (for 
constant m and constant q), several aspects of the resulting metric remain to be investigated. 
In particular, the sources of the electromagnetic field deserve further analysis in this case. 

Since we are mostly interested in the electromagnetic sources, we write here only the 
metric contributions to the Einstein tensor which do not depend on H(t ), 


G t s~ir 

t — Lr r — — Lt q — 


_ r v g 2 ft) 

v a 2 {t)r 4 ' 


(33) 


Note that these are exactly the same as for the charged SD metric (cf. equations (28)). These 


terms should be canceled by the relative components of the electromagnetic stress-energy 
tensor (see next). 
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6.2. Faraday and stress-energy tensors, conserved electric charge and current-density 

Applying the analysis of section [2] to the Thakurta metric (311 we find the same Faraday- 
Maxwell tensor as for the SD case, F tr = —F rt = E(r,t) = qo h(t) /(r 2 a 3 (t)) , where, as 
pointed out above, qo is a constant and h(t) is an arbitrary function of the cosmological 
time t alone, and the only nonzero component of the current-density is J r (r,t) = 
—qo h(t)/(Air a 3 (t)r 2 ). The electromagnetic stress-energy tensor has the following nonzero 
terms 

i]h 2 {t) 


rpt _ rpr _ r?0 

hi t — hi r — -C/0 


= -E* = - 


87 t o 4 (t)r 4 


(34) 


Interestingly, the electromagnetic field and the related stress-energy tensor of the charged 


Thakurta metric (311 are quite the same as for the charged SD geometry (271. 

The components (34) should cancel exactly the associated terms in the Einstein 


tensor (33). Hence, we must choose %h 2 {t) = q 2 (t)a 2 (t). As in the case of the SD metric (see 


section|5f7 the choice of constant q(t) = qo as done in m implies in h(t) = a(t ) and leads to the 
presence of a global electric flow in the radial direction with J r (r,t) = —qoH(t)/(4:ira 2 (t)r 2 ), 
and the total charge of the source varies with time, Q(t) — qoa(t). On the other hand, 
the choice Q = constant (h(t) = 1) has the required physical property that the total electric 
charge of the source is constant with the cosmological time. This corresponds to q(t) = Q/a(t ), 
resulting in a different metric from the solution presented in Ell¬ 
in conclusion, as in the SD case, there are two interesting cases to consider. The first case 
is the one with h{t) = a(t), q(t) = qo = constant, in which the total charge of the source is 
Q(t ) = qoa(t), and there is a nonzero radial electric current J r {r,t). As a matter of fact, this 
is the charged Thakurta solution investigated in (2F, but the presence of the nonzero radial 
electric current was not reported in that study. The second case is the one with h(t) = 1, 
Q(t ) = qo = constant, and q(t) = qo/a(t), what leads two a new metric of Thakurta type, 
which is briefly investigated below. 

To complete this section let us comment about the changes in the energy-momentum 
tensor implied by the replacement of constant m and q by the functions m(t) = M/a(t) and 
q[t) = Q/a{t). The time-time component of the Einstein equations yields G t t — 8 nE t t = 

- + fS) ~ p) ’ where /(r - () = 1 - 2M/a{t)T+ e 2 /"^)-- 2 - The firet 

term has the same form as in the original Thakurta metric, but now with f(r,t) replacing 
f(r) of the original solution. Considering the source as a fluid with heat flow, we find that the 
energy density of the fluid acquires extra contributions which changes the functional form of 
the energy density. The original energy density is positive in all region of the spacetime where 
f(r) > 0, while the new energy density becomes negative for values of R = a(t)r very close to 
the largest zero of the function /(r, t). Similar changes occur in the pressure and heat flow, 
but we do not show such functions here. This negative energy density may be interpreted as 
a kind of phantom matter. On the other hand, we can avoid such a situation by considering, 
for instance, a mixture of fluids instead of a single fluid with heat flow. Other interesting 
way out of this problem is to take different functions for m(t) while keeping q(t) = Q/a(t), 
which is the appropriate electromagnetic source. To simplify the analysis and for the sake of 
comparison with the McVittie metric, we treat here the case m(t) = M/a(t) only. 


6.3. Singularities and horizons of the charged Thakurta metric 

As in the previous cases, for simplicity, we comment only on the undercharged case. 
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For constant m and constant g, the Ricci scalar for the metric (31) may be writte n as 

The 


1Z = 12H 2 (t)/f(r) + 6H(t)/f(r), and the Kretschmann scalar is given in Appendix A 


other curvature scalars present the same singularities as those found from Ricci (A.9) and 


Kretschmann (A.10). 


As seen from the Ricci scalar there are singularities at regions of the spacetime where at 
least one of the following conditions is satisfied: when H{t) diverges, when a(t) vanishes, 
and where f{r) = 0. Thence, besides a possible big-bang singularity when a[t ) —► 0 
(H(t) —> oo), the singularities are defined by the zeros of function f(r) = 1 — 2 m/r + q 2 /r 2 , 
i.e., at r = r± = m ± \]m 2 + q 2 . Of course, the existence of such singularities depends 
upon the asymptotic behavior of the Hubble parameter H(t) for t —> oo. For instance, if 
lim^oo H(t ) =0 and liin^oo H(t) = 0, then the singularity at f(r) = 0 that exists for finite 
t disappears in the infinite future of time t. 

Furthermore, from the Kretschmann scalar (and from other curvature scalars such as 
R^ V R^ V ) it follows that r = 0 is also a singularity. Hence, the singularity is at r = 0 in the 
case f(r) has no real roots, i.e., in the overcharged case, where the central object has electric 
charge q larger than the mass parameter m, to 2 < q 2 . 

The same singularities are present for non-constant m(t) and q(t). The only difference 
is that the singularities defined by the zeros of the function f(r,t) = 1 — 2 m(t)/r + q 2 (t)/r 2 
now depend on time. The curvature scalars for this case are not written here since their 


expressions are cumbersome (see Appendix A). 


As done in the previously analyzed spacetimes, it is useful to define a physical radial 


coordinate by R = r a(t ), so that metric (31) assumes the form 


ds 2 = - f(R) - 


H 2 (t) Rf 

1W 


l2 2 H(t)R , dR 
dt 2 - dtdR + 


f(R) 


f(R) 


R 2 dQ 2 , 


(35) 


where f(R) = 1 — 2 M/R + Q 2 /R 2 and the new parameters M = m(t ) a(t) and Q = q(t) a(t) 
were defined. Again, the search for apparent horizons is performed just in the case 
m(t) = M/a(t ) and q(t) = Q/a(t), with constant M and Q. 

In terms of the new parameters M and Q , for M 2 > Q 2 the curvature singularities 
are located at loci of the spacetime where f(R) = 1 — 2 M/R + Q 2 /R 2 = 0, i.e., R± = 
M ± \JM 2 — Q 2 , where R± = a(t) r± = a(t) (m ± \/m 2 — q 2 ^. In the overcharged case, 
M 2 < Q 2 , there is a curvature singularity at R = 0. 


Since there are singularities in the generalized Thakurta spacetime given by metric (35), it 


is interesting to verify whether there are horizons hiding the singularity to external observers, 


so we look for apparent horizons. In terms of the metric (351, apparent horizons, if exist, are 
given by the roots of the equation H 2 (t) R 2 — f 2 {R) = 0, or 


F T (R,t)=H 2 (t)R 2 - + 

which results in a third-order polynomial equation for R. 


& 

R 2 


= 0 , 


(36) 


The analytical expressions for the three roots of the above polynomial equation (361 for 


R ( R\(t ), R 2 (t), and Rs(t)) are not written explicitly here. The solutions depend explicitly on 
the Hubble parameter H(t). Restricting the analysis to expanding universes H(t) > 0, so that 
H(t) R > 0, the nature of such roots depends on the charge to mass relation. A deeper analysis 
shows that for M 2 > Q 2 one of the roots, Ri(t), say, is smaller than AK = M — \JM‘ 2 — Q 2 
for all times, while the other two are larger than R + = M + \JM 2 — Q 2 . On the other 
hand, for M 2 < Q 2 there is only one real positive root of equation (36), but as mentioned 
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above we report here mostly on the undercharged case (M 2 > Q 2 ). The apparent horizons of 
the extremely charged Thakurta metric (M 2 = Q 2 ) have a similar time behavior as for the 
undercharged case. 

If M 2 > Q 2 there is a curvature singularity at R = R + , with finite t, and the three roots 
may be real and positive. According to the Big-Bang cosmological models, H(t ) is large at 
early times and then, for small times, equation (36) has only one real positive root (R[ (f)), the 
other two roots being complex conjugate to each other. As time goes on, H(t) decreases and 
there is a specific time tf, for which 27 (—1 + 6 MH(tb) — Q 2 H 2 (tb)) = 2 (1 — 6 MH(tb)) 3 ^ 2 
when these roots become real and equal to each other, Ra{tb) = Rs{tb)- Thereafter, for all 
later times, the three horizons develop independently. The situation is depicted in figure [4j 
where we plot R + and the three solutions of equation (36) as a function of the time coordinate 
t for two cases of a(t). At early times, the singularity R = f?+ (t finite) is completely naked, 
and after some intermediate time t = t^ an apparent horizon is formed at R = 
it immediately bifurcates into two branches that evolves independently. 



Figure 4. The evolution of the apparent horizons as a function of the time t in the 
undercharged case of the charged Thakurta spacetime ( |35| >, for two different expansion 
factors, as in the Vaidya (figure |2|) and SD (figure |3| metrics. The plots are for M = 1.0, 
Q = 0.5, a(t) given by equation with k = 0.1 for case (a), and a(t) given by 
equation o with a = 2/3 for case (b). The initial singularity at R = R+, t finite, 
is indicated. For late ti mes , there are three real positive roots of equation FT{R,t) = 0, 
with F t (R, t) given by ( |36| ). 


The solution (i) = R 3 (t), which is satisfied at a particular time, bifurcates and after 
that i ?2 ( t ) and R 3 ( t ) evolve forming a cosmological apparent horizon (or Hubble type surface) 
i? 3 (f), and a black hole type apparent horizon R 2 (t). Asymptotically (for t —>• 00), in 
the case (a) where lim^c*, H(t) = H a > 0, R^it) tends to a de Sitter type cosmological 
horizon, lim^oo R^it) = R^ c = constant in time (i? 3 C is the largest real root of the equation 
1 — 2 M/R + Q 2 /R 2 — H 0 R = 0). In the case (b), where Hindoo H(t ) = 0, Rs(t) tends to the 
Hubble radius, f? 3 (f) ~ 1/t. On the other hand, in both cases (a) and (b), I? 2 (t) tends to a 
black hole horizon. In the case (a) lim^oo I? 2 (t) = f? 2 c = constant (which is the intermediate 
solution of the equation 1 — 2M/R + Q 2 /R 2 — H 0 R = 0), while in case (b) R 2 (t) tends to the 
Reissner-Nordstrom black hole horizon R + = M + \JM 2 — Q 2 . Following [16], see also [18], 
we see that the region (R 2 , t = 00) is a black hole horizon R. Notice also that, in the limit 
t — > 00, for case (b) the Ricci scalar tends to zero and the regions of the spacetime R + and 
R _ are no longer curvature singularities. For this particular case, the true singularity is at 
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R = 0, as it follows from the Kretschmann scalar (A.10). With this we may infer that the 
boundaries of the resulting spacetime are similar to the Reissner-Nordstrom case. 

Lets stress also that if the universe has a phase dominated by the cosmological constant 
(at late times) the Hubble parameter H(t) tends to a constant, and i? 2 (t) tends to a constant 
value larger than R + , i.e., lim^oo i? 2 (t) = Ric = constant, and thus a black hole horizon is 
formed. The cosmological-like horizon Rs(t) tends to a de Sitter horizon (see [16]). In this 
case we may infer that the causal structure of the region R > R + of the charged Thakurta 
type spacetime (351 is similar to the McVittie metric mm- 


7. Further comments and conclusion 

The McVittie metric was proposed in 1933 and have being studied by many authors. Though, 
we can say that a good comprehension of the global structure of the corresponding spacetime 
has been achieved only recently. It is then clear that a lot of effort is still needed to reach a 
good level of understanding on the other metrics that have been written to represent black 
holes in cosmological backgrounds. This is so for the Vaidya metric j2U], and also for the 
Sultana-Dyer [22] and Thakurta [19] metrics. What we mean is that not even the original 
versions of such metrics were investigated in detail, and so any progress in the charged cases 
is of high interest. 

We started our study reviewing the charged metrics given in m , and soon verified that 
the electromagnetic field of those solutions deserved further examination, since in some cases 
the correct sources of the Maxwell electromagnetic field were not given. In particular, in the 
charged versions of the SD and Thakurta metrics the presence of a radial current density 
across the spacetime was not observed. Then, by the sake of comparison, we started by 
investigating this issue in the charged McVittie solution given in [25], and also considered 
a charged Vaidya type metric. In both cases, Maxwell equations were solved by choosing 
the Faraday-Maxwell field in such a way that the four-current density vanishes everywhere 
except at singular points (regions) of the spacetime. The total charge of the central object 
is a constant, and the stress-energy tensor emerging from such a Faraday-Maxwell gauge 
field exactly cancels the corresponding terms in the Einstein tensor produced by the given 
metric. On the other hand, in the charged versions of the SD and Thakurta type metrics, 
the situation is a little more interesting. The Faraday-Maxwell tensor field related to zero 
four-current outside the source contributes to a stress-energy tensor which does not cancel the 
corresponding terms of the Einstein tensor coming from the metrics given in [21] . This gives 
rise to additional energy and pressure terms that have no direct physical interpretation. On 
the other hand, if one chooses the correct Faraday-Maxwell field so to get rid of the undesirable 
energy terms, then there is a global (radial) electric current throughout the spacetime, and 
the electric charge of the source varies with the cosmological time. 

As commented above, further studies are necessary to complete the analysis of the metrics 
considered here. As we have shown, some of them present properties that may represent 
charged black holes in expanding universes. It is worth mentioning once more that the 
geodesic structure of most of those spacetimes were not investigated in the literature, not 
even of the uncharged original metrics. We are now investigating the global structure of the 
charged SD and Thakurta metrics. The complete analysis, including the conformal diagrams, 
will be published in separate papers. 
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Appendix A. Curvature scalars 
a) Charged McVittie metric: 

The charged McVittie metric gets a simpler form by using a Schwarszchild-like radial 
coordinate R given by 


R = a(t) r g(r 7 1 ) = a(t) r ( 1 + 


m 


2 a(t) \ 


4 a(t) \ 


(A.l) 


In terms of this new radial coordinate, and taking constant mass and charge parameters, 
m = M and q = Q, the Ricci and the Kretschmann scalars, 1Z and /C, respectively, for the 
charged McVittie metric (20 1 may be written as 

6 H(t) 


1Z = 12 H\t) + 


JC =48 


VJW) 

M^_Q_ 

R 3 R 4 


+ % + 2W\t) + ^L 
R S U y/m 


3H(t) + 6H 2 {t) - 


2 Q 2 

R 4 


(A.2) 


(A.3) 


where we defined f(R) = 1 — 2M/R + Q 2 /R 2 . These curvature scalars are used in the study 
performed in section [3] 

b) Charged Vaidya type metrics 


The Ricci and the Kretschmann scalars, 7Z and /C, respectively, for the Vaidya metric (22), 
considering m = constant and q = constant, are given by 


1Z = 2 H 2 {t) 6 + 


9(r) 


+ 


2 H{t) (3 


+ 


2 g(r)\ 4 mH(t) 


\ a 2 {t) J a 3 (t)r 2 


a 2 {t), 

1C = 12 (W 4 (f) + 2 H 2 (t)H{t) + H 2 {t )) + ( H A (t) + 2 H 2 (t)H{t) + 2H 2 (t)) 

9 2 (r) 


(A.4) 


+ 


a 2 (t) 


rr2,.\ , o- 4 / + \ 16 H 3 (t)\ 1 6g(r)H(t) , 2 


8 H z (t) + H i (t) - 


8 H(t)q 2 16 mH(t) 

a 4 (t)r 4 
16 H(t) f 6 q 


?(t)7 


a 3 (t)r 2 

1 Qrn n 2 


i[t)r ) 
(H(t)-H 2 (t)) + 


\t)r 


( q 2 — mr ) [1 — a(t)rH(t)\ 


8 H 2 (t) f 7q' 


a e (t)r 2 


20 mq 2 ^ 16m 2 


13mq 2 8 m 


+ 


+ 


\t)i 


7 q 4 12 mq 


3 ' r 2 

2 6m 2 ' 


(A.5) 


2m q 2 

where g(r) = - 77 . These scalars are used in the analysis of the singularities of the 

r r z 

charged Vaidya spacetime, in section [4] 

Now, taking time dependent mass and charge parameters, i.e., with m = m{t) and 
q = q{t ), besides the terms showed by equation (A.41 and (A.5) with m and q replaced 
by rn(t) and q(t), the resulting Ricci and Kretschmann scalar acquire a large number of 
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additional terms. Thus, to avoid cumbersome expressions in the paper, we do not write them 
here. Such extra terms contain negative powers of a{t ), positive powers H(t) and H(t), and 
negative powers of r. These extra terms contain also positive powers of m(t ), g(t), m(t), and 
q{t). 

If we take m(t) = M/a(t) and q(t) = Q/a(t), with constant M and Q, as considered in 
the text, the Ricci scalar becomes 

R = 12iC( f )+3ff«) (2+$j|) + A(1)<?2 
2 M Q 2 


(A.6) 


where g(r,t) = 


a(t)r a 


a 4 (t)r 2 ’ 

-. The Kretschmann scalar results a too long expression to be 


z 2 (t)r 2 

written here, but it has the same kind of terms as in the case of constant m and q. Hence, 
the singularities are the same as in the case with constant m and q. Namely, a big-bang type 
singularity when a(t) = 0, or, equivalently, when H(t) and 77(f) become infinitely large, and 
a singularity at the region (r = 0, t finite). 


c) Charged Sultana-Dyer type metrics 

For constant m and g, the Ricci and the Kretschmann scalar obtained from metric (271 
are, respectively, 

6 (2iC(l) + m) 


n = < 


K = 12 


(1 + g(r)) 2 (2 17 4 (f) + 2 H 2 (t)H(t) + H 2 (t ) 


a(t)r : 

16 g 2 17(f) 


(A.7) 


a 3 (f)r 7 


(2 q 2 — 3mr ) 


1 6H(t)H{t) 
a(t)r 3 
32 H 2 (t) 
a 2 (t)r 6 


( q 2 g ( r ) + Tnr) — 


8 H{t) 
a 2 (t)r 4 


(i + sM) - 


48 mH 3 (t) 
a(t)r 2 


(! + g{r)) 


(g 4 — 3mg 2 r + 3 m 2 r 2 ) + 


a 4 (t)r s 


(7g 4 — 12mg 2 r + 6 m 2 r 2 ) , 


(A.8) 


2 m g 2 
where g(r) = - 

The modifications of these scalars in the case m = M/a(t) and g = Q/a(t ), with con¬ 
stant M and Q , correspond to the presence of new terms with the same functional depen¬ 
dence on a(t), H(t ) , H{t), and r as the above terms. For instance, the Ricci scalar changes 

/ 2 M Q 2 \ ■ ( 5 M Q 2 \ 8 MH(t) 

to TZ = 6 H 2 (t) 2 + 2 r-+ 2 H{t) 3 + —-r--— y~. We do 

w V a(t)r a 2 {t)r 2 J w \ a{t)r a 2 {t)r 2 J a(t)r 2 

not write the resulting expression for the Kretschmann scalar because it is too cumbersome. 

These scalars are used in the analysis of singularities of the charged SD spacetime in section[5j 

d) Charged Thakurta type metrics 

Taking constant parameters m and g, the Ricci and the Kretschmann scalars of the 


charged Thakurta metric (31) may be written respectively as 

617(f) 


12 H 2 (f) 

' f(r) f(r) ’ 


(A.9) 


K = 


+ 


12 

/ 2 ( r ) 


(H(t) + H 2 {t)) 


48 

« 4 (t) 


8q 2 H{t) 12 H 4 (t) 

a 2 (t)r 4 /(r) + / 2 (r) 

16 H 2 (t) \ 8g 4 

a 2 (f) / 2 (r )) a 4 (t) r 8 


(A. 10) 
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These scalars are used in the analysis of section [6] 

As in the SD metric, replacing m and q by function of time m{t) and q does 
note affect the big-bang singularity a(t) = 0, or H(t ) = oo and H(t ) = oo. The 
singularities at the zeros of the function f(r,t) = 1 — 2 m(t)/r + q 2 {t)/r 2 are also 
present, but now depend on the cosmic time t. We do not write here the ex¬ 
pressions for the Ricci and Kretsclimann scalars for the general because they are too 
large. In the particular case where m(t) = M/a(t) and q(t) = Q/a(t), the Ricci 
scalar becomes 1Z = 12 H 2 (t)/f(r,t) — 2 H(t) (6 f(r,t) + M/a(t)r — Q 2 /a 2 (t)r 2 ) / f 2 (r,t) + 

2 H 2 [t) A(M/a{t)r — Q 2 /a?(t)r 2 ) 2 — f(r,t) (6M/a(t)r —5 Q 2 /a 2 (t)r 2 ) /f 3 (r,t), and we 
avoid to write the huge expression for the Kretschmann scalar. 
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